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Abstract 

We report improved measurements of B decays with an n meson in the final state using 357 fb _1 
of data collected by the Belle detector at the KEKB e + e~ collider. We observe the decays B^ — > 
rjir^ and B^ — > nK^; the measured branching fractions are B(B ± — > vtt^ 1 ) = (3.9 ± 0.5(stat) ± 
0.2(sys)) x 1CT 6 and B{B ± -> nK^ = (2.2 ± 0.4(stat) ± O.l(sys)) x 10 -6 . Their corresponding 
CP-violating asymmetries are measured to be —0.10 ± O.ll(stat) ± 0.02(sys) for nir^ and —0.55 ± 
0.19(stat)^QQ3(sys) for r/K^. No significant signal is found for B° — ► nK° decays; the upper limit 
on the branching fraction at the 90% confidence level is B(B° — > rjK°) < 1.9 x 10~ 6 . 

PACS numbers: 13.25.Hw, 12.15.Hh, 11.30.Er 
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Charmless B decays provide a rich sample to understand B decay dynamics and to 
search for CP violation. It has been suggested [l( that direct CP violation could be large 
in B + — ► t]tx + j2j and B + — > r]K + decays due to penguin-tree interference. In our earlier 
measurements |3[, no CP-violating asymmetry (Acp) was seen for the t]ti + mode. Belle's 
central value of Aqp for r]K + was large in magnitude and negative, but statistics were 
limited even for the branching fraction measurement. In the most recent update from the 
BaBar collaboration the previous negative Acp values with ~ 2er significance for both 
r]K + and r]Ti + modes have moved to within ~ la deviation from zero. Because experimental 
uncertainties on CP asymmetries and branching fractions are still large, it is necessary to 
analyze larger data sample. 

In this paper, we report improved measurements of branching fractions and partial rate 
asymmetries for B — > rjh decays, where h is a K or ir meson. The partial rate asymmetry 
for charged B decays is defined to be: 

N(B~ -> •qhr) - N(B + -> nh + ) 
CP ~ N{B- -> nh-) + N{B+ r,h+) ' ^ 

where N(B~) is the yield for the B~ — > r]h~ decay and N(B + ) denotes that of the charge 
conjugate mode. The data sample consists of 386 million BB pairs (357 fb" 1 ) collected 
with the Belle detector at the KEKB e + e~ asymmetric-energy (3.5 on 8 GeV) collider [(| 
operating at the T(4S) resonance. 

The Belle detector is a large-solid-angle magnetic spectrometer that consists of a silicon 
vertex detector (SVD), a 50-layer central drift chamber (CDC), an array of aerogel threshold 
Cerenkov counters (ACC), a barrel-like arrangement of time-of-flight scintillation counters 
(TOF), and an electromagnetic calorimeter (ECL) comprised of CsI(Tl) crystals located 
inside a superconducting solenoid coil that provides a 1.5 T magnetic field. An iron flux- 
return located outside of the coil is instrumented to detect K\ mesons and to identify muons 
(KLM). The detector is described in detail elsewhere 0. In August 2003, the three-layer 
SVD was replaced by a four-layer radiation tolerant device. The data sample used in this 
analysis consists of 140 fb _1 of data with the old SVD (Set I) and 217 fb _1 with the new 
one (Set II). 

The event selection and candidate reconstruction are similar to that documented in our 
previous paper Here, we only give a brief description of event selection. Two r\ decay 
channels are considered in this analysis: rj — > 77 (f/ 77 ) and rj — > 7r + 7r~7r° (773^). We require 
photons from rj and ir° candidates to have laboratory energies above 50 MeV. In the r/ 77 
reconstruction, the energy asymmetry, defined as the absolute value of the energy difference 
in the laboratory frame between the two photons divided by their energy sum, must be 
less than 0.9. Neither photon is allowed to pair with any other photon with at least 100 
MeV energy to form a 7r° candidate. Candidate 7r° mesons are selected by requiring the 
two-photon invariant mass to be in the mass window between 115 MeV/c 2 and 152 MeV/c 2 . 
The momentum vector of each photon is then readjusted to constrain the mass of the photon 
pair to the nominal ir° mass. 

Candidate r/ 37r mesons are reconstructed by combining a 7r° with at least 250 MeV/c 
laboratory momentum with a pair of oppositely charged tracks that originate from the 
interaction point (IP). We make the following requirements on the invariant mass of the 77 
candidates in both data sets: 516 MeV/c 2 < M 77 < 569 MeV/c 2 for r/ 77 and 539 MeV/c 2 < 
M% n < 556 MeV/c 2 for 773^. After the selection of each candidate, the r\ mass constraint is 
implemented by readjusting the momentum vectors of the daughter particles. 
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Charged tracks are required to come from the IP. Charged kaons and pions, which are 
combined with rj mesons to form B candidates, are identified using a K(tt) likelihood LK{L n ) 
obtained by combining information from the CDC (dE/dx), the TOF and the ACC. Dis- 
crimination between kaons and pions is achieved through a requirement on the likelihood 
ratio Lk/ (L w + L^). Charged tracks with likelihood ratios greater than 0.6 are regarded as 
kaons, and less than 0.4 as pions. Furthermore, charged tracks that are positively identified 
as electrons or muons are rejected. The K/tt identification efficiencies and misidentification 
rates are determined from a sample of D* + — > D°ir + ,D — > K~tt + decays. In reference 
[H we reported that in both data sets n + has a higher efficiency than 7r~ while the kaon 
identification efficiency is higher for K~ than K + . The efficiency difference introduces a 
bias in Acp which needs to be corrected. Kg candidates are reconstructed from pairs of 
oppositely-charged tracks with an invariant mass (M„) between 480 and 516 MeV/c 2 . Each 
candidate must have a displaced vertex with a flight direction consistent with that of a Kg 
originating from the IP. 

Candidate B mesons are identified using the beam constrained mass, M^ c = 



y-E'beam — -fflj an d the energy difference, AE = Eb — -Ebcam, where -Ebcam is the run- 

dependent beam energy in the T(4S') rest frame and is determined from B — > D^n events, 
and Pb and Eb are the momentum and energy of the B candidate in the T(4S') rest frame. 
The resolutions on M^ c and AE are about 3 MeV/c 2 and 20-30 MeV, respectively. Events 
with Mbc > 5.2 GeV/c 2 and \ AE\ < 0.3 GeV are selected for the analysis. 

The dominant background comes from the e + e~ — > qq continuum, where q = u, d, s or c. 
To distinguish signal from the jet-like continuum background, event shape variables and B 
flavor tagging information are employed. We combine information of correlated shape vari- 
ables into a Fisher discriminant [9J] and compute the likelihood as a product of probabilities 
of this discriminant and cos#b, where 9b is the angle between the B flight direction and the 
beam direction in the T(4S') rest frame. A likelihood ratio, 1Z = C S /(C S + C qq ), is formed 
from signal (C s ) and background {C qq ) likelihoods, obtained using events from the signal 
Monte Carlo (MC) and from data with M^ c < 5.26 GeV/c 2 , respectively. Additional back- 
ground discrimination is provided by B flavor tagging. An event that contains a lepton (high 
quality tagging) is more likely to be a BB event so a looser 1Z requirement can be applied. 
We divide the data into six sub-samples based on the quality of flavor tagging (see de- 
tails in [3]). Continuum suppression is achieved by applying a mode dependent requirement 

on 1Z for events in each sub-sample in Set I and Set II according to A^r xp / J Ns xp + Ng^ p , 
where N^ xp is the expected signal from MC and A^| p denotes the number of background 
events estimated from data. 

From MC all other backgrounds are found to be negligible except for the r]K + <-> f]ir + 
reflection, due to K + <-> tt + misidentification, and the feed-down from charmless B decays, 
predominantly B — > rjK* (892) and B — > ?7p(770). We include these two components in the 
fit used to extract the signal. 

The signal yields and partial rate asymmetries are obtained using an extended unbinned 
maximum-likelihood (ML) fit with input variables M^ c and AE. The likelihood is defined 
as: 




C = exp 



n(£NjPj) and 



(2) 



-[l-qi-Acp^PjiM^AEi), 



(3) 
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where % is the identifier of the i-th event, P(M\ )C , AE) is the two-dimensional probabil- 
ity density function (PDF) in Mb c and AE, q indicates the B meson flavor, B + (q = +1) 
or B~(q = —1), Nj is the number of events for the category j, which corresponds to ei- 
ther signal, qq continuum, a reflection due to K-n misidentification, or background from 
other charmless B decays. For the neutral B mode, Vj in the equation above is simply 
Pj(Mb ci , AEi) and there is no reflection component. 

Since the efficiency of particle identification is slightly different for positively and neg- 
atively charged particles, the raw asymmetry defined in Eq. 1 must be corrected. This 
efficiency difference results in an Acp bias of —0.005 (0.005) for tj-k (rjK). The bias is 
subtracted from the value obtained for the raw asymmetry. 

The signal peak positions and resolutions in Mb c and AE are adjusted according to the 
data-MC differences using large control samples of B — > Dir and D° — > K + ix~n° /7r°7r° 
decays. The continuum background in AE is described by a first or second order poly- 
nomial while the Mb c distribution is parameterized by an ARGUS function, f{x) = 
xy/1 — x 2 exp [— £(1 — x 2 )}, where x is M^ c divided by half of the total center of mass 



energy Thus the continuum PDF is the product of an ARGUS function and a polyno- 
mial, where £ and the coefficients of the polynomial are free parameters. The PDFs of the 
reflection and charmless B backgrounds are modelled as two-dimensional smooth functions, 
obtained using large MC samples. 

All the signal and background yields and their partial rate asymmetries are free param- 
eters in the fit except for the reflection components, where the Acp and the normalizations 
are fixed to expectations based on the B + — > t]K + and B + — > t]tx + partial rate asymmetries 
and branching fractions, as well as K + <-> tt + fake rates. The reflection yield and Aqp are 
first input with the assumed values and are then recalculated according to our measured 
results. 

Table U shows the measured branching fractions for each decay mode as well as other 
quantities associated with the measurements. The efficiency for each mode is determined 
using MC simulation and corrected for the discrepancy between data and MC using the 
control samples. Other than the particle identification performance discrepancy reported 
earlier 0], our MC slightly overestimates the detection efficiency of soft n°s in Set II, which 
results in a 4.3% correction for the t] 3n mode. The combined branching fraction of two data 
sets are computed as the sum of the yield divided by its efficiency in each set divided by the 
number of B mesons. The combined branching fraction of the two r\ decay modes is obtained 
from the weighted average assuming the errors are Gaussian. Systematic uncertainties in 
the fit due to the uncertainties in the signal PDFs are estimated by performing the fit after 
varying their peak positions and resolutions by one standard deviation. In B ± — > r]ir ± , the 
reflection yields in Set I (Set II) are estimated to be 3.4 (5.5) events for the ?? 77 mode and 
1.1 (2.0) for r] 3n . And in B — > f]K ± , the reflection yields are 3.2 (6.5) for r^ 77 and 1.3 (2.2) 
for 77 37r . The reflection yields and their Aqp values are varied by one standard deviation 
in the fit to obtain the corresponding systematic uncertainties. The quadratic sum of the 
deviations from the central value gives the systematic uncertainty in the fit, which ranges 
from 3% to 6%. For each systematic check, the statistical significance is taken as the square 
root of the difference between the value of —2 ln£ for zero signal yield and the best-fit value. 
We regard the smallest value as our significance including the systematic uncertainty. The 
number of B + B~ and B°B° pairs are assumed to be equal. 

The performance of the 1Z requirement is studied by checking the data-MC efficiency 
ratio using the B + — > D°tt + control sample. The obtained systematic error is 2.0-3.3%. 
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TABLE I: Detection efficiency (e) including sub-decay branching fraction, yield, significance (Sig.), 
measured branching fraction (B), the 90% C.L. upper limit (UL) and Acp for the B — > rjh decays. 
The first errors in columns 3, 5 and 7 are statistical and the second errors are systematic. 



Mode e(%) Yield Sig. £(l(r 6 ) UL(10" 6 ) A 



CP 



J B ± ^r ? vr ± 10.9 3.9 ±0.5 ±0.2 -0.10 ± 0.11 ± 0.02 

?7 77 7r ± 8.7 3.9tgj±0.2 0.08 ± 0.14 ± 0.02 

Set I 9.12 71.5±J|;1 ± 2.0 7.0 5.2+^ ± 0.3 0.14 ± 0.18 ± 0.02 

Set II 9.34 68.3±i£;|±|;5 5.1 3.1 ±0.7 ±0.2 0.00 ± 0.22±g;g| 

7)3^ 6.7 3.8±g;| db 0.3 -0.54 ± 0.22 ± 0.02 

-6.6 j_ n c q q q 0+1-3 _i_ n q nmxn qe+0.02 

0.03 



Set I 3.26 16.1Ij$±0.8 3.3 3.2+^ ± 0.3 -0.10 ±0.38] 

Set II 3.53 33.7£gi2i° 5.8 ^Lo ± °- 3 -O^et^-aw 
B^^-qK* 7.1 2.2 ±0.4 ±0.1 -0.55 ± 0. 19^;^ 



^ 5.8 2.0 ± 0.5j£? -0.7llg^;SS 

Set I 8.35 27.0^° °1 2 2 3 1 3.3 2.1±g;f ± 0.2 — 0.44±g;|f ± 0.03 

Set II 8.47 39.0±il;?±|:| 4.5 2.0t££ ± 0.2 -0.94^ ± 0.07 

7f± 4 12 s+°' 8+0 - 2 _n 28+°- 29 + - 04 

'/37r-"- z -°-0.6-0.1 u - z °-0.31-0.03 

Set I 3.03 7.3±t:t±o:2 22 1-6±S;| ± 0.2 — 0.66±?;fJ±g;ig 

Set II 3.19 21.2i££tgJ 3.5 2.8±J;g ± 0.2 -0.20±^J| ± 0.03 

B°^ V K° 2.1 0.9 ±0.6 ±0.1 < 1.9 

?? 77 /f° 1.0 0.6±g;J ±0.1 < 2.2 

Set I 3.06 -1.6i|3+o ; 5 _ — o.3±g;? ± 0.1 < 1.7 

Set II 3.17 9.0t£§ ± 0.7 1.4 1.2±^ ± 0.1 < 3.1 

?? 37r K° 1.8 l-5±i:g ±0.1 < 3.9 

Set I 0.90 3.5±||±g| 1.5 2.3±f;| ± 0.3 < 6.0 

Set II 1.00 2.4±|;f±g;| 1.0 l.OtJjj ± 0.1 < 3.7 



The systematic errors on the charged track reconstruction are estimated to be around 1% 
per track using partially reconstructed D* events. The tt° and ?7 77 reconstruction efficiency 
is verified by comparing the tt° decay angular distribution with the MC prediction, and by 
measuring the ratio of the branching fractions for the two D decay channels: D° — > K + tt~ 
and D° —>■ K + ir~ir . We assign a 3.5% (4.0 %) error for the ir° and 7/ 77 reconstruction in 
Set I (Set II). The Kg reconstruction is verified by comparing the ratio of D + — > K s ir + and 
D + — > K ~7r + 7r + yields. The resulting K s detection systematic errors for Set I and Set II 
are 4.4% and 4.0%, respectively. The uncertainty in the number of BB events is 1%. The 
final systematic error is obtained by first summing all correlated errors linearly and then 
quadratically summing the uncorrelated errors. 

Figure^ shows the M^ c and AE projections after requiring events to satisfy —0.10 GeV 
< AE < 0.08 GeV and Mb c > 5.27 GeV/c 2 , respectively. No significant signal is observed 
for the B° — > t]K° meson mode. Significant signals are observed for charged B decays; the 
corresponding M^ c and AE projections for B + and B~ samples are shown in Figs. El and 
01 The Aqp results for the two rj decay modes in two data sets are combined assuming that 
the errors are Gaussian. No significant asymmetry is observed for the mode while an 
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FIG. 1: M hc and AE projections for (a,b) B ± r ? vr ± , (c,d) B ± -> ijK^, and (e,f) B° r]K° 
decays with the rjyy and 773^ modes combined. Open histograms are data, solid curves are the fit 
functions, dashed lines show the continuum contributions and shaded histograms are the feed-down 
component from charmless B decays. The small contributions around M^ c = 5.28 GeV/c 2 and 
AE = ±0.05 GeV in (a)-(d) are the reflection backgrounds from i? ± — > r\K^ and i? ± — > tjtt^. 

excess of B + over B~ decays is seen in the r/K^ mode, giving a large negative Aqp value. 

In summary, we have observed — > tjtt^ and B^ — > r]K ± decays; their branching 
fractions are measured to be (3.9±0.5±0.2) x 10~ 6 and (2.2±0.4±0.1) x 10~ 6 , respectively. 
The central values of our measurements in both modes are 1.4<r smaller than the results [4| 
reported by the BaBar collaboration. Our Aqp measurement for B ± — > 7777 ,— 0.10±0.11± 
0.02, is consistent with no asymmetry. However, the partial rate asymmetry for B ± — > rjK^ 
remains large in magnitude and negative, Acp = —0.55 ± 0.19l~o;o3, which is 2.9<r from zero. 
Larger data samples are needed to verify this large CP asymmetry. Finally, no significant 
B° — > rjK° signal is found and we assign an upper limit at the 90% confidence level of 
B(B° ^r]K°) < 1.9 x 10" 6 . 

We thank the KEKB group for the excellent operation of the accelerator, the KEK 
Cryogenics group for the efficient operation of the solenoid, and the KEK computer group 
and the National Institute of Informatics for valuable computing and Super-SINET network 
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FIG. 2: Afbc an d AE projections for (left) B~ — > r]ir~ and (right) B + — > 7771" + with the ry 77 and 773^ 
modes combined. Open histograms are data, solid curves are the fit functions, dashed lines show 
the continuum contributions and shaded histograms are the contributions from charmless B decays. 
The small contributions near = 5.28 GeV/c 2 and AE = —0.05 GeV are the backgrounds from 
misidentified B^ — > rjK^ (reflections). 
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FIG. 3: Mfr c and AE projections for (left) B~ — > r]K~ and (right) B + — > r]K + with the r]~ n and 
7737T modes combined. Open histograms are data, solid curves are the fit functions, dashed lines 
show the continuum contributions and shaded histograms are the contributions from charmless B 
decays. The small contributions near M^ c = 5.28 GeV/c 2 and AE = 0.05 GeV are the backgrounds 



from misidentified B^ — ► (reflections). 
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